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(54) Title: METHOD AND APPARATUS FOR SEGMENTATION OF AN OBJECT 




(57) Abstract: The invention relates to a method and an apparatus for segmentation of an object in a 2D or 3D image data set by 
extracting a path along the object. In order to obtain a higher accuracy and reliability, the method according to the invention comprises 
the steps of: a) selecting a start point (5) of the path (2) as first active point, b) adapting an adaptable model (40) to the object (1) in 
a first active section around the start point (5), c) finding the next point (58) of the selected region by the steps of: cl) copying the 
geometric model parameters of the adapted model (40) of the active point to a plurality of neighboring points (51-58) of the active 
point (5), c2) orienting a model (32) around each of said plurality of neighboring points (51-58) using said copied model parameters 
and, searching the closest object points (21-26) around each of said neighboring points (51-58) based on said model, c3) adapting 
the models to the found object points (21-26) for each neighboring point (51-58), c4) selecting the neighboring point (58) for which 
the adapted model (48) fits best to the object (1) as next point of the selected region and as next active point, and d) repeating step c) 
until an end point of the path or a predetermined number of iterations is reached. The invention may preferably be applied for vessel 
segmentation to create maximum intensity projections with overlapping/occluding vessels removed. 
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Method and apparatus for segmentation of an object 



The invention relates to a method and an apparatus for segmentation of an 
object in a 2D or 3D image data set by extracting a path along the object in a selected region. 
Further, the invention relates to a computer program product. 

Magnetic resonance angiography (MRA) images provide important 
5 information for the diagnosis of vascular disease, such as arterial stenosis and aneurysm. The 
recent development of MR blood-pool contrast agents which have extended intravascular 
half-life allows the acquisition of high resolution, high contrast 3D images of the vascular 
system. 

The visualization of the vessel pathways is crucial to allow quick and reliable 

10 assessment of any potential problems. The most common visualization method is to construct 
a maximum intensity projection (MIP). However, with blood-pool contrast agents, the longer 
scan times necessary to achieve higher resolution require imaging during the steady state of 
contrast agent diffusion. Therefore, both arteries and veins are enhanced, and diagnostically 
important information (typically the arteries, where stenosis occurs) may be fully or partially 

1 5 occluded in the MIP. 

Several approaches have been made for selecting voxels belonging to vessel 
regions. A wide spread approach for vessel enhancement is to use multi-scale orientation 
selective filters, based on eigen-analysis of the Hessian matrix, as e.g. known from "Model- 
based quantitation of 3D magnetic resonance angiographic images", A. Frangi et al., IEEE 

20 Transactions on medical imaging, Vol. 18, No. 10, October 1999. Therein, linear vessel 

segments are modeled with a central vessel axis curve coupled to a vessel wall surface. The 
path is initialized using the shortest path across an image iso-surface. However, a significant 
limitation of filter-based approaches for vessel selection is that not all those voxels required 
to define a vessel structure fulfil the filter criteria, in particular those pixels near structural 

25 bifurcations. On the other hand, pixels not belonging to vessels may also be selected with 
filtering based approaches, for example in regions between two nearby vessel structures. 

Also for the visualization of other objects, particularly anatomical structures 
like bones or airways in medical image data sets and also in 2D image data sets like images 
of the vessel tree, a method for segmentation of an object is often applied. It is thus an object 
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of the present invention to provide a method and an apparatus for accurately and 
automatically segmenting an object in an image data set. 

This object is achieved by a method for segmentation as claimed in claim 1 
comprising the steps of: 
5 a) selecting a start point of the path as first active point, 

b) adapting an adaptable model to the object in a first active section around the 
start point, 

c) finding the next point of the selected region by the steps of: 

cl) copying the geometric model parameters of the adapted model of the 
10 active point to a plurality of neighboring points of the active point, 

c2) orienting a model around each of said plurality of neighboring points 
using said copied model parameters and searching the closest object points around each of 
said neighboring points based on said model, 

c3) adapting the models to the found object points for each neighboring 

15 point, 

c4) selecting the neighboring point for which the adapted model fits best 
to the object as next point of the selected region and as next active point, and 

d) repeating step c) until an end point of the path or a predetermined number of 
iterations is reached. 

20 According to the invention the user selects a start point as first action point in 

a first active region of the selected region in which the path shall be extracted. Preferably, 
also an end point is selected. Thus, the object of interest is identified which is then 
automatically segmented so that it can be suppressed from an image if required. For said 
segmentation prioritized region growing is used wherein iteratively voxels/pixels are added 

25 to a selected region based on a model adaptation and a selection for which point the model 
fits best to the object. The method according to the invention thus includes an iterative 
algorithm for finding the points of the path along the object during which the model, which 
has been adapted around a previous point of the path, is first copied to all neighboring points. 
Thereafter, for each neighboring point the models are adapted to the object by finding the 

30 closest object points. Based on the selection which model of which neighboring point fits 
best to the object the next point of the path from said neighboring points is selected. 

During said prioritized region growing geometric model parameters, e.g. the 
radius, or cross-section parameters of the model, are propagated to neighboring points not yet 
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having parameter estimates. In the first step, this means all the neighbors, but in subsequent 
steps, some of the neighbors will already have parameter estimates. 

By the method according to the invention the boundary of the object can be 
accurately identified. Further, since particularly in vascular images venous and arterial 
5 pathways are often close together, it is possible, via the use of an appropriate geometric 
model, to discriminate between very closely separated structures, so that only anatomically 
connected pathways are selected. Further, the method is able to detect objects across a range 
of scales, which is important since the width of objects like vessels or airways can vary 
significantly. 

10 Preferred embodiments of the invention are included in the dependent claims. 

The object is also achieved by an apparatus for segmentation as claimed in claim 12 
comprising: 

a) start point selection means for selecting a start point of the path as first active 

point, 

15 b) first adaptation means for adapting an adaptable model to the object in a first 

active section around the start point, 

c) path extracting means for finding the next point of the selected region by the 
steps of: 

cl) copy means for copying the geometric model parameters of the 
20 adapted model of the active point to a plurality of neighboring points of the active point, 

c2) search means for orienting a model around each of said plurality of 
neighboring points using said copied model parameters and for searching the closest object 
points around each of said neighboring points based on said model, 

c3) second adaptation means for adapting the models to the found object 
25 points for each neighboring point, 

c4) neighbor selection means for selecting the neighboring point for 
which the adapted model fits best to the object as next point of the selected region and as next 
active point, and 

d) control means for repeating step c) until an end point of the path or a 
30 predetermined number of iterations is reached. 

The invention relates further to an apparatus for acquiring and processing 
medical image data, in particular a magnetic resonance apparatus, a computer tomography 
apparatus, an x-ray apparatus or an ultrasound apparatus, comprising means for acquiring 
medical image data and means for processing said image data including an apparatus for 
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segmentation according to claim 13. Still further, the invention relates to a computer program 
product comprising computer program means for causing a computer to perform the steps of 
the method as claimed in claim 1 when said computer program product is run on a computer. 
It shall be understood that these apparatuses and said computer program product can be 
5 developed further and can have similar embodiments as included in the dependent claims 
referring to claim 1 . 

For the selection which model fits best to the object, i.e. which of the 
neighboring points is the next point of the path, a distance measure is used based on the 
distances between the object and the model. Therefore, the distances along the normals from 

10 the surface of the model to the found object points are formed, and the distances along said 
normals are processed to form a distance measure for each model. This distance measure is 
combined with a gradient measure (sum of gradient across cylinder surface), and the model 
having the minimum combined measure is selected. The corresponding neighboring point 
will be selected as next point to be added to the selected region. 

15 In order to find object points in the image data set that are closest to the 

models around the new neighboring points gradients of image values in the image data set are 
used. Particularly in medical images objects like vessels or airways show a high gradient at 
their borders compared to the gradient of the grey values of surrounding tissue. 

In general, the invention can also be used to segment a surface. Therefore, the 

20 surface is extracted by using planes along the surface as adaptable models and by adapting 
the planes to the surface. 

A preferred application of the invention lies, as explained above, in the field of 
segmentation of objects in medical images. Preferably, tubular objects like vessels, bones and 
airways can be segmented in 3D medical image data sets. The start point is then selected 

25 within the tubular object, and the path within the object can be extracted, the path following 
mainly the centerline of the tubular object. 

In general, any kind of three-dimensional model can be used when segmenting 
an object in a 3D image data set. When applying the invention for segmentation of tubular 
objects, cylinder models have been shown useful the cross section of said cylinders being 

30 either circular, elliptical or even more close to the cross section of the tubular object to be 
segmented. 

When using such cylinder models, in the first step of region growing, when 
parameters of the cylinder model are copied to all neighboring points of the previous point of 
the path, preferably the orientation and the radius of the cylinder model is copied to each 
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neighboring point so that identical cylinder models are located around each neighboring 
point. After finding the closest object points to such cylinder models, the cylinder models are 
then adapted to the object by adapting the orientation and the radius of the cylinder models 
according to the detected object points. 
5 According to the invention a cylinder model based boundary response is 

directly incorporated as a priority function during region growing leading to an increased 
accuracy of the segmentation. After completely extracting the path it is then possible to refine 
the models by combining all the models to create a single deformable model for the selected 
vessel segment, ant then to refine its surface, particularly by using a mesh on the surface and 
10 adapting the mesh to the object using a known method. 

Preferably, the path is recovered from all selected points of the selected region 
by following said points in the order in which they were added to the selected region. 



1 5 The invention will now be explained in more detail with reference to the 

drawings in which: 

Fig. 1 shows a flow-chart of the method according to the invention, 
Fig. 2 illustrates the main steps of the method, 
Fig. 3 shows cylinder model adaptation, 

20 Fig. 4 shows a first step of prioritized region growing, 

Fig. 5 shows another step of prioritized region growing, 
Fig. 6 illustrates a step of finding the next point of the path, 
Fig. 7 shows region growing according to the invention, 
Fig. 8 shows the adaptation of cylinder models along a path, 

25 Fig. 9 illustrates the use of the method for segmentation of a surface, and 

Fig. 10 shows a medical imaging apparatus according to the invention. 



The flow-chart shown in Fig. 1 illustrates the main steps of the method for 
30 segmentation according to the invention. In a first step SI the user selects a start point for a 
path along the object to be segmented thereby also indicating the object that is to be 
segmented. Said start point is used as first active point in the selected region. In a second step 
S2 an adaptable model, e.g. an adaptable geometric primitive cylinder model with adaptable 
radius and orientation, is arranged around the start point and is adapted to the object. For 
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adaptation object points, often also referred to as feature points, are searched for around the 
start point so that the model fits as good as possible to the local intensity distribution. 

In order to find the next point of the selected region an iterative algorithm 
comprising several steps S3-S11 is performed. 
5 In a recursion comprising steps S3 to S 10 a model shall be adapted to each 

neighboring point of the first active point. The first step S3 of said recursion checks if the 
recursion is finished, i.e. if a model has been adapted for all neighboring points of the first 
active point. If this is the case a jump is made to step S9 where from all said neighboring 
points the next point of the selected region is selected. If not yet all neighboring points are 

10 done in step S4 the next neighboring point is selected to which a model shall be adapted. 

Thereafter, in step S5, it is checked if said selected neighboring point has already geometric 
parameters for arranging a model around it or not. If the neighboring point has already 
parameters a jump is made to step S7 while in the other case the parameters of the adapted 
model of the first active section, i.e. of the section around the start point (the first active 

15 point) are copied to the present neighboring point. When using a cylinder said model is 
characterized by its length, its radius and its axial orientation. Around a neighboring point 
thus an identical model is arranged using the parameters copied from the active point. 

In step S7 the closest object points are searched for the present neighboring 
point. Based on the found object points a model arranged around the neighboring point in 

20 step S6 is now adapted to the found object points, i.e. the size and orientation of said model is 
adapted such that it fits as good as possible to the object points found for the present 
neighboring point (step S8). In order to repeat these steps for all neighboring points of the 
present active point the algorithm jumps back to step S3. 

If for all neighboring points a model has been adapted such that each 

25 neighboring point has an individual model, the neighboring point is selected as next point of 
the selection region based on said individual models which fits best to the object (step S9). 
The selected neighboring point is then used as next active point, and the section around said 
neighboring point will be used as next active section if - after a negative result of the 
decision in step S10 - the algorithm continues with step S3 in order to find the successive 

30 next point in the selected region. 

If an end point has been reached or if a predetermined number of iterations of 
said recursive algorithm has been made the query of step S10 gives a positive result so that 
the method for segmentation continues with a final step SI 1 in which all points in the 
selected region can be combined in the order they have been found thus giving the path along 
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the object. All adapted models of said points can then be combined and refined using a mesh 
on the surface of said models, particularly a triangular mesh which can be more closely 
refined to the object using a known method. The obtained segmentation data may then be 
used to suppress the object from a maximum intensity projection obtained from a 3D image 
5 data set if the object occludes any other information included in the 3D image data set but not 
visible in the maximum intensity projection. 

The method according to the invention can be implemented such that after 
selection of an object of interest the object is segmented and suppressed automatically. 
Compared to known methods using multi-scale filtering the method according to the 

10 invention achieves an increased accuracy during segmentation so that a final image has a 
higher image quality, and the probability that the object is incorrectly segmented is 
considerably reduced. 

Particular steps of the method according to the invention will now be 
illustrated by use of a particular example showing vessel segmentation in maximum intensity 

1 5 projections for magnetic resonance angiography using cylinder models. Since arteries and 

veins are often fully or partially occluded in a maximum intensity projection image it is often 
required to select either of them for segmentation in order to suppress it from the projection 
image. Typically, stenosis occurs in the arteries so that veins should be segmented and 
suppressed. The main steps will be hereinafter explained. 

20 Fig. 2 schematically illustrates the main steps of the method for segmentation 

according to the invention. Fig. 2 A shows how adaptable cylinder models 3 are adapted 
within the vessel 1 according to boundaries of the vessel 1 (the object). After a series of such 
adaptations, the points with the best fitting models may be extracted to obtain the path 2, 
together with associated radius and orientation estimates. This is illustrated in Fig. 2B. Fig. 

25 2C shows a mesh 4 constructed using the information from the cylinder models 3 which may 
be deformed to refine the final segmentation. This is done after the path 2 has been 
completely extracted and all cylinder models 3 have been found and adapted to the object 1, 
in order to even more adapt the cylinder models 3 to the object. As a result, segmentation 
data for the object 1 are obtained. 

30 Fig. 3 illustrates the adaptation of a cylinder model 3. Assuming that there is a 

start point 5 selected within the tubular object 1, i.e. within the vessel 1, an initial cylinder 
model 30 arranged around said start point 5 (first active point) in the start section is adapted 
to the boundary 17 of said vessel 1. The length, radius and orientation of the adapted cylinder 
model 40 is found by searching object points 11-16 starting from the initial cylinder model 
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30 by searching along the normals 6 on the surface of said cylinder model 30. Object points 
1 1-16 are particularly characterized by a significant increase in the gradients of the image 
values, i.e. the grey values of an image data set, compared to the gradients of the tissue inside 
or outside the boundaries 17 of the vessel 1. When such object points 1 1-16 are found, the 
5 initial cylinder model 30 is adapted to fit as good as possible with such object points 11-16 
resulting in the adapted cylinder model 40. Such adaptation step is thereafter also used during 
the recursive algorithm in order to find the next points in the selected region. 

In Figs. 4-6 the steps for finding the next point in the selected region (steps 
S3-S9 in Fig. 1) shall be illustrated in 2D for illustration only, while in reality the cylinder is 

10 fitted in 3D. Starting from a start point 5 for which an adapted cylinder model 40 has already 
been found having a parameter set P 4 o one of the neighboring points 51-58 has to be selected 
as next point of the path. At first the parameters P 4 o of the adapted cylinder 40 are copied to 
all neighboring points 51-58, and around each neighboring point 51-58 an identical cylinder 
is arranged, i.e. the parameter sets P 3 i to P38 for cylinder models 31 to 38 (not shown) 

15 arranged around each of said neighboring points 51-58 are identical to the parameter set P40 
of the cylinder model 40 arranged around the start point 5. 

In a next step shown in Fig. 5 the parameter sets P31 to P38 are updated into 
parameter sets P41 to P 4 8 by transforming the original identical cylinder models 31-38 into 
adapted cylinder models 41-48 being adapted to the object. It should be noted that before 

20 adaptation of the cylinder models 3 1-38 it is required for each neighboring point 51-58 to 

find the closest object points so that the cylinder models 31-38 can be adapted to such object 
points. 

Fig. 6 illustrates the steps using the example shown in Fig. 3. It should be 
noted that the dimensions and distances are not to scale in order to improve distinctiveness. 

25 In reality the distances between neighboring points are much smaller. 

The start point 5 together with its adapted cylinder model 50 as well as all the 
neighboring points 51 to 58 are shown. As explained above, around each neighboring point 
51-58 a cylinder model 31-38 identical to the cylinder model 40 is arranged. For simplicity's 
sake only one cylinder model 32 around neighboring point 52 is shown. Next, for each of 

30 said cylinder models 31-38 the closest object points of the surface 17 of the object 1 in 
directions orthogonal to the surfaces of the cylinder models 31-38 are searched. Using 
cylinder model 32 as an example it can be seen that on one side (left side in the drawing) 
object points 1 1, 12, 13 along the surface 17 of the object 1 can be found while on the other 
side (right side in the drawing) no object points can be found. While this will be similar for 
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the cylinder models 31, 33, 35, 36, 37 around the neighboring points 51, 52, 53, 55, 56, 57, 
located outside the object 1, it will be different for the cylinder models 34, 38 arranged 
around the neighboring points 54 and 58, respectively, which are located inside the object 1. 
For those cylinder models 34, 38 in all directions perpendicular to the cylinder axis object 
5 points, e.g. objects points 21-26 for cylinder model 38 around neighboring point 58, can be 
found. 

Next, the original cylinder models 31-38 are adapted to the object using the 
found object points for each cylinder model, i.e. the size and orientation of said cylinder 
models is changed so as to fit each cylinder model as good as possible to the found object 

10 points. This is illustrated by way of the exemplary cylinder models 42, 47 and 48. Cylinder 
model 42 is fit to the object points 11, 12, 13 by changing its radius and, eventually, its 
orientation, cylinder model 47 is adapted to the object 1 to fit as good as possible to the 
object points 24, 25, 26, and the cylinder model 48 is fit to the object points 21-26 by not 
only amending their radius but also the orientation around its neighboring point. As can 

15 immediately be seen cylinder model 48 fits best to the object 1 - provided that the path to be 
extracted shall go into the direction of neighboring point 58 and not into the direction of 
neighboring point 54. 

As a measure for the decision which of the adapted models 41-48 fits best with 
the object a distance measure is preferably used. For each of said adapted cylinder models the 

20 distances between the adapted cylinder model and the corresponding found object points in 
directions perpendicular to the cylinder axes are measured. Such distances, e.g. distances d 2 i 
to d26 for the cylinder model 48, are then used to determine a distance measure for said 
cylinder model, e.g. by calculating the mean square root of said distances. Since for cylinder 
models 41-43 and 55-57 not enough distances can be determined since not in all directions 

25 perpendicular to the cylinder axes object points can be found, a distance measure for such 

cylinder models cannot be calculated or is assumed to be very high so that the corresponding 
neighboring points 51-53 and 55-57 are automatically excluded from being one point along 
the path. As next point along the path the neighboring point 58 is selected which has the 
minimum distance measure. Thereafter, the same recursive algorithm is used to find the next 

30 point along the path from all neighboring points of the next active point 58 until an end point 
is reached or until a predetermined number of iterations of said recursive algorithm has been 
gone through. 

After object point detection, cylinder parameters can be updated as follows. A 
new axis orientation can be determined as the mean orientation over all vectors between 
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object points at opposite ends of the cylinder for the same radial orientation. The updated 
radius of the cylinder can then be calculated as the mean perpendicular distance of the object 
points to the updated axis. 

In the particular application of the invention for segmenting a vessel as criteria 
5 for extracting points of the path inside the vessel a vesselness response for the adapted 

cylinder models can be defined using the residual distances between the object points and the 
adapted cylinder surface as well as individual feature strength leading to a speed function in 
which certain parameters control the sensitivity of surface evolution to the respective terms. 

The described front propagation approach is well suited to selecting vessel 

10 structures where a simple local structure is repeated to form a complex pattern at larger 

scales. Initial orientation parameters are estimated at the start point via an exhaustive search, 
then propagated as voxels (in 3D image data sets) moved into a border set, and updated 
whenever the time is computed. A path will then be constructed between the start point and 
an end point by following the order in which points are added. 

15 Fig. 7 shows how several points 5a-5f are extracted to subsequently form the 

path along the object. Fig. 8 shows an extracted path. Around the points 5, 5a, 5b, . . ., 5f the 
adapted cylinder models 60, 71-76 and, for all point except the start point 5, the non-adapted 
cylinder models 61-66 the parameters of which are copied from the cylinder model 60-64 of 
the respective previous point 5, 5d, are shown. It should be noted again that the object is 

20 typically much larger than the voxel spacing and that there are many more intermediate steps 
than those shown. However, these intermediate steps cannot be easily illustrated, since the 
individual cylinders of neighboring voxels overlap significantly. 

According to the invention the centerline extracted via a prioritized region 
growing is used to reconstruct the vessel volume. If a cylinder-based speed function is used, 

25 orientation and radius estimates are available directly, otherwise cylinder models may be 
oriented along the path adapted to retrieve these estimates. However, a vessel's cross 
sectional profile often deviates from circular. Visualization applications require accurate 
detection of the vessel wall, in order to avoid residual regions appearing in maximum 
intensity projections. Therefore, preferably a deformable model is constructed using the 

30 centerline and radius estimates, which can be adapted to refine the segmentation. 

The vessel boundary estimate can, for instance, be represented using a 
triangulated mesh, which is adapted according to image features while also imposing shape 
based constraints on the deformation. The mesh is formed by placing a set of disks at discrete 
steps along the path. The radius and orientation of each disk are determined as the average of 
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values of the available estimates for all points within the neighborhood of the disk. A set of 
points around the circumference of each disk is then used to define the triangulated mesh. 

The chosen mesh adaptation method maintains the underlying triangle 
structure and uses the initializing configuration as a guiding shape model to avoid excessive 
5 deformation. Adaptation is an iterative procedure consisting of a surface detection, similar to 
the object point detection described above, followed by minimization of an energy function. 
The energy is composed of an external, image-related energy and an internal, shape-related 
term wherein the relative influence of each term is weighted. The external energy attracts the 
mesh towards the surface points. The internal energy is defined with respect to changes of the 

10 difference vectors between neighboring mesh vertices, penalizing large deviations from the 
initial shape. Energy minimization uses the conjugate gradients method. 

The invention can also be applied for segmentation of a surface 8. This is 
illustrated in Fig. 9. Therein planes E0-E3, shown schematically in 2D cross-section for 
simplicity, are used as models to be adapted to the surface 8. Starting from a start point 80 for 

15 which a plane E0 is adapted to the surface 8 the parameters of the plane E0 are copied to all 
neighboring points 81, 82, 83 through which identical planes (not shown) are arranged. Based 
on such planes the closest object points, i.e. points on the surface 8, are searched in directions 
perpendicular to the planes. Using the found object points the initial planes through the 
neighboring points 81, 82, 83 are then adapted resulting in the shown adapted planes El, E2, 

20 E3. Based on these adapted planes El, E2, E3 the distances between the planes and the 
surface 8 are then determined in directions perpendicular to the planes, and a distance 
measure for each neighboring point 81, 82, 83 is calculated. As can immediately be seen in 
Fig. 9 the neighboring point 81 will have the minimum distance measure, thus leading to the 
selection of neighboring point 81 as next point in the selected region and as next start point 

25 for finding the next point in the region. 

A medical imaging apparatus including an apparatus for segmentation 
according to the invention is shown in Fig. 10. Therein image acquisition means 91 are used 
for acquiring medical images of an object of interest 90, which maybe a patient from which 
angiographic image data shall be obtained. The image acquisition means 91 may include any 

30 kind of medical imaging modality like magnetic resonance, computer tomography , x-ray or 
ultrasound. The acquired image data are then inputted into image processing means 92 
including a segmentation unit 94 and a post-processing unit 95. The obtained image data may 
then displayed on a display 93. 
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The method for segmentation as explained above will be performed by the 
segmentation unit 94. Said unit includes a start point selection unit 941 for selecting a start 
point of the path and a first adaptation unit 942 for adapting adaptable model to the object in 
a first active section around the start point. A path extracting unit 943 is thereafter used to 
5 extract all the points along the path. Said unit 943 includes a copy unit 945 for copying the 
geometric model parameters of the adapted model of the start section to a plurality of 
neighboring points of the start point, particularly those points which have no current 
parameters, and orienting a model around each of said neighboring points using said copied 
parameters, a search unit 946 for searching the closest object points around each neighboring 

10 point based on said models, a second adaptation unit 947 for adapting the models to the 

found object points for each neighboring point and a neighbor selection unit 948 for selecting 
the neighboring point for which the model fits best to the object as next point in the selected 
region and as next active point. In addition, control means 944 are provided for controlling 
the method for segmentation until an end point or a number of iterations is reached. If the 

1 5 extraction of the path is finished, in the post-processing unit 95 the obtained adapted models 
can be refined as explained above using a mesh, and the segmented object can be removed 
from the initial image data set, if required, in order to make objects visible in the image data 
set, e.g. in a maximum intensity projection, that had been occluded previously by the 
segmented object. 

20 The invention has been explained by using vessel segmentation as a particular 

application. However, the invention is not limited to such application. The invention may 
also be applied to segment other objects like airways, bones or, preferable and more general, 
tubular objects. 

The invention may further be applied for segmentation, in 2D or 3D image 
25 data sets, or other objects, via the use of other geometrical primitives, e.g. plane. In the case 
of the cylinder model, using the invention has been shown to provide an increased accuracy 
leading finally to a higher image quality and reliability. 

The invention uses prioritized region growing, where selection of points to 
include into the region is achieved on the basis of a "goodness of fit 11 response for a 
30 geometrical shape primitive (or model). Each point in the region has an associated shape 
model. Model parameters are propagated from the selected region to other boundary points 
(not selected, but bordering the selected region), as initial estimates for parameters at these 
point, which are then updated via feature search and model fitting using the new feature 
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points. After model fitting during the parameter update step, the "goodness of fit" measure 
can be obtained by comparing the fitted model and the extracted feature points. 
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1 . Method for segmentation of an object in a 2D or 3D image data set by 

extracting a path (2) along the object (1) in a selected region, the method comprising the 
steps of: 

a) selecting a start point (5) of the path (2) as first active point, 

5 b) adapting an adaptable model (40) to the object (1) in a first active section 

around the start point (5), 

c) finding the next point (58) of the selected region by the steps of: 

cl) copying the geometric model parameters of the adapted model (40) of 
the active point to a plurality of neighboring points (51-58) of the active point (5), 
10 c2) orienting a model (32) around each of said plurality of neighboring 

points (51-58) using said copied model parameters and searching the closest object points 
(21-26) around each of said neighboring points (51-58) based on said model, 

c3) adapting the models to the found object points (21-26) for each 
neighboring point (51-58), 
1 5 c4) selecting the neighboring point (58) for which the adapted model (48) 

fits best to the object (1) as next point of the selected region and as next active point, and 

d) repeating step c) until an end point of the path or a predetermined number of 
iterations is reached. 

20 2. Method according to claim 1, characterized in that the selection in step c4) 

which model fits best to the object (1) is based on a distance measure of the distances (d 2 i- 
d26) between the object (1) and the model (48). 

3. Method according to claim 1, characterized in that the search of the object 

25 points (21-26) in step c2) is based on the gradients of image values in direction of the model's 
surface normal in the image data set. 
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4. Method according to claim 1, characterized in that the closest object points 

(21-26) are searched in directions perpendicular to the orientation copied to each neighboring 
point (51-58). 

5 5. Method according to claim 1 , characterized in that the method is used for 

segmentation of a surface (8) and that for extracting the path along said surface planes (E0- 
E3) are used as adaptable models. 

6. Method according to claim 1, characterized in that the method is used for 

10 segmentation of tubular objects in a 3D image data set, in particular tubular medical objects 
like vessels, bones, airways in a 3D medical image data set, and that for extracting the path 
inside such a tubular object a 3D adaptable model is used. 

7. Method according to claim 6, characterized in that a cylinder model (40) is 
15 used as adaptable model. 

8. Method according to claim 7, characterized in that in step cl) the orientation 
and the radius of the cylinder model (40) is copied to each neighboring point (51-58) and that 
a cylinder model (32) having said copied orientation and radius is arranged around each 

20 neighboring point (51-58). 

9. Method according to claim 1, characterized in that the method is used for 
segmentation of vessels in magnetic resonance angiography, in particular for removing the 
segmented vessels from maximum intensity projections. 

25 

10. Method according to claim 1, characterized in that the adapted models along 
the extracted path are combined and/or refined using a mesh on the surface of said models for 
refinement of the segmentation. 

30 11. Method according to claim 1, characterized in that the path is recovered from 

all selected points of the selected region by following said points in the order in which they 
were added to the selected region. 
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12. Apparatus for segmentation of an object in a 2D or 3D image data set by 

extracting a path (2) along the object (1) in a selected region, the apparatus comprising the 
steps of: 

a) start point selection means (941) for selecting a start point (5) of the path (2) 
5 as first active point, 

b) first adaptation means (942) for adapting an adaptable model (40) to the object 
(1) in a first active section around the start point (5), 

c) path extracting means (943) for finding the next point (58) of the selected 
region by the steps of: 

10 cl) copy means (945) for copying the geometric model parameters of the 

adapted model (40) of the active point to a plurality of neighboring points (51-58) of the 
active point (5), 

c2) search means (946) for orienting a model (32) around each of said 
plurality of neighboring points (51-58) using said copied model parameters and for searching 
15 the closest object points (21-26) around each of said neighboring points (51-58) based on said 
model, 

c3) second adaptation means (947) for adapting the models to the found 
object points (21-26) for each neighboring point (51-58), 

c4) neighbor selection means (948) for selecting the neighboring point 
20 (58) for which the adapted model (48) fits best to the object (1) as next point of the selected 
region and as next active point, and 

d) control means (944) for repeating step c) until an end point of the path or a 
predetermined number of iterations is reached. 

25 13. Apparatus for acquiring and processing medical image data, in particular 

magnetic resonance apparatus, computer tomography apparatus, x-ray apparatus or 
ultrasound apparatus, comprising means (91) for acquiring medical image data and means 
(92) for processing said image data including an apparatus (94) for segmentation according to 
claim 12. 

30 

14. Computer program product comprising computer program means for causing a 

computer to perform the steps of the method as claimed in claim 1 when said computer 
program product is run on a computer. 
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